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Abstract 

We analyze the transport properties of contacts between unconventional superconductor and nor- 
mal diffusive metal in the framework of the extended circuit theory. We obtain a general boundary 
condition for the Kcldysh-Nambu Green's functions at the interface that is valid for arbitrary 
transparencies of the interface. This allows us to investigate the voltage-dependent conductance 
(conductance spectrum) of a diffusive normal metal (DN) / unconventional singlet superconductor 
junction in both ballistic and diffusive cases. For d-wave superconductor, we calculate conductance 
spectra numerically for different orientations of the junctions, resistances, Thouless energies in DN, 
and transparencies of the interface. We demonstrate that conductance spectra exhibit a variety 
of features including a F-shaped gap-like structure, zero bias conductance peak (ZBCP) and zero 
bias conductance dip (ZBCD). We show that two distinct mechanisms: (i) coherent Andreev re- 
flection (CAR) in DN and (ii) formation of midgap Andreev bound state (MABS) at the interface 
of d-wave superconductors, are responsible for ZBCP, their relative importance being dependent 
on the angle a between the interface normal and the crystal axis of d-wave superconductors. For 
a = 0, the ZBCP is due to CAR in the junctions of low transparency with small Thouless energies, 
this is similar to the case of diffusive normal metal / insulator /s-wave superconductor junctions. 
With increase of a from zero to 7r/4, the MABS contribution to ZBCP becomes more prominent 
and the effect of CAR is gradually suppressed. Such complex spectral features shall be observable 
in conductance spectra of realistic high-Tc junctions at very low temperature. 

PACS numbers: PACS numbers: 74.20.Rp, 74.50.+r, 74.70.Kn 
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I. INTRODUCTION 



The low energy transport in mesoscopic superconducting systems is governed by Andreev 
reflection a unique process specific to electron scattering at normal metal/superconductor 
interfaces. The phase coherence between incoming electrons and Andreev reflected holes per- 
sists at a mesoscopic length scale in the diffusive normal metal, which enhances interference 
effects on the probability of Andreev reflection P| . The coherence plays an important role at 
sufficiently low temperatures and voltages when the energy broadening due to either voltage 
or temperature becomes of the order of the Thouless energy Exh of the mesoscopic structure. 
As a result, the conductance spectra of mesoscopic junctions may be significantly modified 
by these interference effects. A remarkable experimental manifestation of the electron-hole 
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phase coherence is the observation of the zero bias conduct ance_pea^ 
normal metal (N) / superconductor (S) tunneling junctions 

Various theoretical models of charge transport in diffusive i unctions extend the clean 
limit theories dev elop ed bv Blonder, Tinkham and Klapwijk Q] (BTK) and Zaitsev. |l^ 



In Refs. 



QHEiE 




2l\ the scattering matrix approach was used. On the other 



hand, the quasiclassical Green's function method in nonequilibrium superconductivity [22] 
is much more powerful and convenient for the actual calculations of conductance for the 
arbitrary bias- voltages Using the Kuprianov and Lukichev (KL) boundary condition 



2J| for a diffusive SIN interface, Volkov, Zaitsev and Klapwijk (VZK) have obtained the 



conductance spectra with ZBCP, origin of which was attributed to coherent Andreev re- 
flection (CAR) which induces the proximity effect in diffusive metal 23|. Several authors 
studied the charge transport in mesoscopic junctions combining this boundary condition 
with Usadel 
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equations that describe superconducting correlations in diffusive metal 
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The modified boundary conditions were studied by several authors j35|, |36[ . Important 
progress was achieved by one of the authors who developed the so-called "circuit the- 

ory" for matrix currents that allows to formulate boundary conditions for Usadel-like equa- 
tions in the case of arbitrary transparencies. By using this generalized boundary condition, 
three of the authors have evaluated the conductance in DN/S junctions and demonstrated 
how various may be the conductance spectra. 38 1. 

Another mechanism of ZBCP plays the role in ballistic unconventional superconductor 
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junctions. The conductance peak in this situation arises from the formation of midgap 
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42| |. The experimental observation 



Andreev bound states (MABS) at the interface 

of the ZBCP has been reported for various unconventional superconductors of anisotropic 

57l |. A basic 



pairing symmetry 
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theory of ballistic transport in the presence of MABS has been formulated in Refs 



42|. 



Stimulated by this theory, extensive studies of MABS in unconventional superconductor 



junctions have been per 
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'ormed during the 
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ast decade: in the case of broken time reversal 



65l |. in triplet superconductor junctions j6 



er effect |75|, 



in quasi-one dimensional organic superconductors 
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79|, MABS in ferro mag net Junctions l80 



l86l . 1871 . |88| , influence of MA BS on Jo s ephson ef fect 



other related problems 
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However, an interesting 



question remained: how two mechanisms of ZBCP (due to CAR and due to MARBS) work 
together, this being relevant for diffusive normal metal / unconventional superconductor 
(DN/US) junctions. 

To solve this problem, three of the present authors have recently extended the c ircui t 



theory to the systems that contain unconventional singlet superconductor junctions [10 
Application of this theory for DN/rf-wave superconductor (DN/d) junctions has shown that 
the formation of MABS strongly competes with the proximity effect that is an essential 
ingredient for CAR in DN. The MABS induces the unconventional channels where quasi- 
particles are resonantly transmitted through the interface. The overall contribution of these 
channels to the proximity effect is however suppressed by the isotropization, i.e., the an gular 
averaging over momentum directions of injected quasiparticles. However, the reference (l08l | 
does not contain the necessary technical details of the matrix current derivation and presents 
the results only for low voltage limit. To compare with experiment, one has to evaluate the 
conductance spectrum in wide range of bias voltage. 

In this paper, we present a detailed derivation of the matrix current in (DN/US) junctions. 
Although the relation obtained is valid for both singlet and triplet superconductor junctions, 
we focus on the case of singlet superconductor. We present detailed numerical calculations of 
the conductance spectra of DN/US junctions for d-wave superconductors. We investigate the 
dependence of the spectra on various parameters: the height of the barrier at the interface, 
resistance Rd in DN, the Thouless energy Eth in DN and the angle between the normal to 
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the interface and the crystal axis of rf-wave superconductor (a). We normahze the voltage- 
dependent conductance as{eV) by its value in the normal state, ctat, so that axieV) = 
as{eV)/aN- 

Our main results are as follows: 

1. The ZBCP is frequently seen in the shape of aT{eV). For a ^ 0, the ZBCP is robust not 
depending on the diffusive resistance Rd- For a = 0, ZBCP is due to the CAR. 

2. The appearance of ZBCP is different for MABS and CAR mechanisms. The first mech- 
anism may lead to arbitrarily large (Jt{0). The second mechanism can not provide o"t(0) 
exceeding unity. While for the first mechanism the width of the ZBCP is determined by the 
transparency of the junction, it is determined by Thouless energy for the second one. These 
two mechanisms comp ete s ince the proximity effect and the MABS in singlet junctions are 
generally incompatible [lOS^ . 

3. In the extreme case a = ir/A the proximity effect and the CAR are absent. The aT{eV) is 
then given by a simple Ohm's lawiarieV) = {Rb + Rd) / {RRa=o + Rd) Rb being the resistance 
of the interface. 

4. For a = 0, when MABS are absent for Rd = 0, the ZBCP of axieV) is attributed to the 

CAR alone. When the transparency of the junction is sufficiently low, axieV) for | eV \< Ag 

is enhanced with the increase of Rd due to the enhancement of the proximity effect (Aq is the 

maximum amplitude of the pair potential). The ZBCP becomes prominent for Eth ^ Aq 

and Rd/Rb < 1 [Rd 7^ 0). In this case, the ZBCP turns into a zero bias conductance dip 

I I 

(ZBCD) with a further increase of Rd/Rb ,34] . 

5. We have shown that the conductance spectrum can vary substantially depending on the 
parameters. The results obtained are important to analyze the actual experimental data on 
conductance spectra of high Tc cuprate junctions since in this case the diffusive scattering 
in the normal metal is of special relevance. 

The structure of the paper is as follows. We formulate the model in use in section 2. 
We also present there the detailed derivation of the matrix current and end up with the 
expression for the normalized conductance. We focus on rf-wave superconductor junctions 
in section 3 and evaluate aT{eV) and the measure of the proximity effect ioi various cases. 
We summarize the results in section 4. 
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II. FORMULATION 



In this section we introduce the model and the formahsm. We consider a junction consist- 
ing of normal and superconducting reservoirs connected by a quasi-one-dimensional diffusive 
conductor (DN) with a length L much larger than the mean free path. The interface between 
the DN conductor and the US (unconventional superconductor) electrode has a resistance 
Rfj while the DN/N interface has zero resistance. The positions of the DN/N interface and 
the DN/S interface are denoted as x = —L and x = 0, respectively. According to the circuit 
theory js^], the constriction area {—Li < x < Li) between DN and US is considered as 
composed of the diffusive isotropization zone {—Li < x < —L2), the left side ballistic zone 
(— L2 < X < 0), the right side balhstic zone (0 < x < Li) and the scattering zone (x = 0). 
The scattering zone is modeled as an insulating delta-function barrier with the transparency 
T„ = 4 cos^ 0/(4 cos^ + Z"^), where Z is a dimensionless constant, is the injection angle 
measured from the interface normal to the junction and n is the channel index. We assume 
that the sizes of the ballistic and scattering zones along x axis is much shorter than the 
superconducting coherence length. 

Here, we express insulating barrier as a delta function model H6{x), where Z is given 
hj Z = 2mH/{h'^kp) with Fermi momentum kp and effective mass m. In order to clarify 
charge transport in DN/US junctions, we must obtain Keldysh-Nambu Green's function, 
which has indices of transport channels and the direction of motion along x axis taking into 
account the proper boundary conditions. For this purpose it is necessary to extend a general 
heory of boundary condition which covers the crossover from ballistic to diffusive cases 
37 1 formulated for conventional junctions in the framework of the circuit theory {2^5 Is? ]. 
However, the circuit theory cannot be directly applied to unconventional superconductors 
since it requires the isotropization. To avoid this difficulty we restrict the discussion to a 
conventional model of smooth interface by assuming momentum conservation in the plane 
of the interface. 

In this section, we will show how to derive the matrix current in DN/US junctions. Then 
we will derive the retarded and Keldysh components of the matrix current. Finally, we will 
show how to calculate conductance of the junctions. 
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A. Calculation of the matrix current 



To derive the relation between the matrix current and Green's functions, we make use 
of the method proposed in ^| . The method puts the older ideas to the framework of 
Landauer-Biittiker scattering formalism. One expresses the matrix current in a constriction 
in terms of one-dimensional Green's functions gn,cr;n',cr'{£',x,x'), where n,n' and a, a' = ±1 
denote the indices of transport channels and the direction of motion along x axis, respec- 
tively. The "check" represents the Keldysh-Nambu structure. These Green's functions have 
to be expressed in terms of the transfer matrix that incorporates all information about the 
scattering, and asymptotic Green's functions presenting boundary conditions deep in each 
side of the constriction. Here, we restrict the discussion to a conventional model of smooth 
interface, assuming momentum conservation in the plane of the interface. Within the model, 
the channel number eventually numbers possible values of this in-plane momentum and the 
transfer matrix becomes block-diagonal in the channel index. Following the treatment de- 
veloped in Ref. 0], we thus solve Green's functions gn,a;n',a'{£,x,x') separately for each 
channel. gn,a;n',a'{£,x,x') can be expressed as 

gn,a;n',a'{£,X,x') = ^ exp{iaPnX - id' PnX')G'^''' {x , x') (l) 
o-,(t'=±1 

The function G'^"' [x, x') are varying smoothly at the scale of 1/p and obey the following 
semi classical equation. 

i^^^n^ + H{x))G''/{x,x') = (2) 
with effective Hamiltonian H[x) where H{x) and Ho{x) are given by 

H{x) = Ho{x) - timp{x), (3) 



Ho{x) = sf, + A(x), = I ° I . (4) 



In the above, T,imp{x) is a self-energy due to the impurity scattering and Sjmp(x) 7^ is 
satisfied only for a: < 0. The self-energy originating from the superconducting pair potential 
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A(x) = for a; < while A (a;) for a; > is given by 

H^)J^''],kJ ° ^"1. (5) 

G^'^' (x, x') has discontinuity at x = x' 

G"/ {x + 0, x) - G"/ (x - 0, x) = -zia5^^/ / I I ■ (6) 
To get rid of the discontinuity, we define (^^""'(x, x') as follows 

2iGl"\x, x') = [gf{x, x') + asign(x - x')]/ \vn\ ■ (7) 

We denote ^^"^'(0+, 0+) = g2 and (7^'^'(0_, 0_) = gi, where g2 and gi satisfy 

^2 = M^^iM, 

using a transfer matrix M. 

The derivation of the matrix current 

/ = ^Tr„,^[E^^i] = ^Tr„,^[S^^2]. (8) 
is given in Appendix, where it is shown that /„o can be represented as follows [lOS^ . 

InO = 2[G'i, Bn] 

with 

5„ = (-Ti„[Gi, HZ'] + HZ'H+ - TIG,HZ'H+Gi)-\TUI - Hz') + T^^G,Hz'H+) (9) 



This is the very equation, which was firstly derived as eq. (2) of Ref. |l08l |. It should be 
remarked that this formula of the matrix current is very general since it is available both for 
singlet and triplet superconductors. For low transparency limit, i.e., T„ << 1, Tin << 1, 
Ino can be approximated to be 

Ino = YK\^-H-),Gi] (10) 

Eq. (fTUj) can be regarded as an extended version of Kuprianov and Lukichev's boundary 
condition 2^ for unconventional superconductor junctions. On the other hand, when G2+ = 
G2- is satisfied as in the conventional superconductor, since 

lim5_^o - D-'[Tin{2Gi - [Hz\ 61]+) + Hz'H+ + TIG^Hz'H+Gi] 



= {H+ + T^nGi)-\-H+ + Ti„Gi) 

and 

lim^__,oGiD-iG'i[Ti„(2G'i - [Hz\Gr]+) +TIHz'H^ + G^Hz'H+Gr] 
= (Gi + TinH+Y^{—TinH+ + Gl) 
are satisfied, then the resulting J„o and / are reduced to be 

4o = 2Ti„(l + T2„ + Ti„[i^+,Gi]+)-i[i7+,G'i] (11) 



n 



which is identical to eq. (36) of Ref. 
Appendix. 



In the above, the definition of D is given in the 



B. Calculation of the retarded part of the matrix current 

In order to calculate the retarded part of the matrix current, we denote Keldysh-Nambu 
Green's function Gi, G'2±, 



G. = ' : , G,.= , (12) 





where the Keldysh component ki^2± is given by ki(2±) = -Ri(2±)/i(2)(0) - /i(2)(0)Ai(2±) 
with the retarded component -Ri,2± and the advanced component Ai^2± using distribution 
function /i(2)(0). In the above, R2± is expressed by 

R2± = i9±T3 + f±f2) 

with g± = ej — A^, /± = A-t/ a/ — e^, and A2± = —t^rI^t^, where e denotes the 
quasiparticle energy measured from the Fermi energy. /2(0) = /os'(O) = tanh[e/(2kBT)] in 
thermal equilibrium with temperature T. Here, we put the electrical potential zero in the 
US-electrode. We also denote H+, H^, B„, I as follows. 
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H^=\ ^ / , iJ-= , (13) 
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Hereafter, in the present paper, we focus on the singlet superconductor junction case without 
broken time reversal symmetry states (BTRSS). For triplet case or singlet one with BTRSS, 



the situation becomes much more complex 



109l | and we will discuss in forthcoming paper 



in detail. In singlet superconductors, we can choose Ri = cos 6*0 T3 + sin 6*0 to satisfy the 
boundary condition at the interface. After some algebra, we can obtain Br as follows 

Br = — Ti„[l + T'}^ + Tin{RiRp^ + Rp^Ri)]^'^[TinRi + Rp^] 

where we have used the relation 

RiR^Rp + R^ RpRi = 
The resulting Ir can be written as 



f 2e2^ -2T„[cosgo(/+ + /_) - singo(^+ + g^)] 

^ (2- TO(l + (?+^?- + /+/-)+ T„[cos^o(^+ + ^-)+sin^o(/+ + ^ ' 

This is one of the central results of this paper 



C. Calculation of the Keldysh part of the matrix current 

Next, we focus on the Keldysh component. We define h 

^>^=jJ2niKf3]. (15) 

n 

After straightforward calculations, h is given by 

h = jYl Trace{f3(i?ifiA' + R\Bk) 

n 

-[fsiRWj, + BrR, + + R\BR)]foNiO) - + rI){Br + 5i)]/3jv(0)}, 

with Ki = ^1/1(0) - /i(0)ii /i(0) = /oAr(O) + /37v(0)f3. Since Br and ^1 is proportional 
to the linear combination of 7^2 and f^, the second term which is proportional to /oAr(O) 
disappears. It is necessary to obtain Bk which is given by 

Bk = Dr'Nk - D-^'DkD-a'Na, (16) 
10 



with 

' Dr Dk 



D = -Ti,[Gi, HZ'] + HZ'H+ - TIG,Hz'H+G,, D 
where Nk and Na is the Keldysh and advanced part of N given by 
N = -Ti^Hz' + T^^GiHz'H+ N = 



Da 



Nr Nk 
Na 



We can express Nk and Dk as hnear combination of distribution functions /o5(0), /oAr(O), 
and /3iv(0) as follows, 

Nk = 6-1/05(0) + C'2/oiv(0) + (73/3iv(0), 

by 2 X 2 matrix d {i = 1..6). Taking account of the fact that Dr'Ci, Dr'C2, D^'CiD^^NA 
and D^C^D'^Na can be expressed by the linear combination of T2 and fa, while D^C^ 
and D]^CqD^ Na are proportional to the linear combination of i and fi, we can express /{, 
as follows, 

2 

h^\Y. 'I^ace{(^i + R\)[BKEh - {Br + ^],)]/3iv(0)}. 



with 



n 

Bke = D-r'[C^-CqD-/NaI 
Cq = Ti„[-(^if3 - f3ii)i-^ + ^-'(^ifa - fsii) 



-T^n{Rih - fsi^i-^ipii - TinRiR-'Rp{Rih - hAi)]. (17) 

Since following equations are satisfied, 

^A^^A = -TsBI^Ts, Am(p) = -f3Rl(j,f3, 
^R^i-TlnRm + TinRlR^Rp) = Br, 

Ijj is given by 

h=^ 5^[-(^i + R\)Br{Ri + Ri)B^R - {Ri + Ri){BR + b],) 
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+Ti„(i?i + Rl)Dn\Ri + Rl){Rl)-'Bi 



+rf„(^i + R\)Dn\Ri + R\)iRl)-'RlR^,Bl]fMO)- (18) 
After a simple manipulation, we can show 

BniRm + TinRlRpRm) ~ TinRm Rp = = cLj^TinRmRp^ , (19) 

with 

, _ (1 + Tl){l + g+g_ + /+/_) + 2Ti4cos^o(g+ + g^) + singo(/+ + /_)] 
^ 1 + ^+5- + /+/- 

Then the resulting is given by 

h^j E[-(^i + R\)MRi + - (^1 + R\){Br + sjj) 

n 

-Tl{R, + ^I)^^^(^i + ^I)(^l,)-^]/3iv(0). (20) 
Since Br is given as Br = h2T2 + h^r^, with 

-Ti„[Ti„ sin ^o(l + 9+9^ + /+/-) + /+ + /-] 



62 = 



(1 + r2J(l + g+g_ + /+/_) + 2ri„[cos^o(^+ + 9-) + sin^o(/+ + /-)] ' 



, ^ -Tin[Tin COS ^0(1 + 9+9- + /+/_) +9+ + 9-] /2ll 

' (l + T2J(l+^+(?_ + /+/_) + 2ri„[cos^o(^+ + ^-) + sin^o(/+ + /-)]' ^ ^ 
final expression of /(, is given by following equation. 

" V 2 I (2-7;)(l + ^+g_ + /+/_)+7;[coseo(^+ + ^-)+sin^o(/+ + /-)] ? 

(22) 

Co = r„(l+ I cos^o r + I sin^o T) 
x[| g+ + g_ p + I /+ + /_ p + I ! + /+/_ +^+^_ p + I Ug.-g+f. \'] 
+2(2 - 7;)Real{(l + g^g*, + /;/!)[(cos^o + cos^*)(^+ + g_) + (sin^o + sin^*)(/+ + /_)]} 
+4rjmag(cos 60 sin ^S)Imag[(/+ + f.){gl + gl)]. 
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This is a very general expression which is available for any singlet superconductor without 
BTRSS. This expression is also one of the central results of this paper. For isotropic limit 
where /+ = /_ and Qj^ = is satisfied, we obtain 

' h 4"2| (2-T„)+T„[^/+cos^o + /+sin^o] r 
Ai = (1+ I cos^^o I' + I sin^^o Hd I' + I /+ P +1) 



+4Imag[/+5(+]Imag[cos 6^0 sin 9*], (23) 



which is identical to eq. (11) of Ref. 



On the other hand for the ballistic limit, where 



6*0 = is satisfied, we can reproduce the generalized BTK formula 



h 2-^ M ^ _ nr.r 12 ^^^> 



T„(l + T„ 


r+ 


' +{Tn - 1) 








1 


+ (T„ - i)r+r_ 


2 



h 

n 

with 



where following equations are satisfied 

9^ ^9- i + r+r_ /+ + /_ 2(r+ + r. 



1 + ^7+^7- + /+/- i-r+r_ 1 + + i-r+r_ 

D. Calculation of the conductance 

In the following, we apply the quasiclassical Keldysh formalism for calculation of the 
conductance. The spatial dependence of 4 x 4 Green's function in DN Gj^six) which is 
expressed in the matrix form as 



Rn{x) Kn{x] 
An{x) 



should be determined. The Keldysh component Kn{x) is given by Kn{x) = RN{x)fi{x) — 
fi{x)AN{x) with retarded component Rn{x), advanced component Am{x) using distribution 
function fi{x). We put the electrical potential zero in the S-electrode. In this case the spatial 
dependence of Gn{x) in DN is determined by the static Usadel equation [2^, 
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with the diffusion constant D in DN, where H is given by 

^0 



H 



Ho 



with Hq — £73. 

The boundary condition for Gat (a;) at the DN/S interface is given by, 

^[G„(.)^l,„o_ = < / > ■ (26) 

The average over the various angles of injected particles at the interface is defined as 

/Tr/2 pn/2 
# cos 0/(0)/ / d<j)T{<j)) COS <j) (27) 
■7r/2 J-Tr/2 

with I {(f)) — I and T{(f)) — Tn- The resistance of the interface Rb is given by 

h 2 



#2"(0) cos(;6 



t/2 

Gn{—L) coincides with that in the normal state. The electric current is expressed using 
Gn{x) as 



I el 



^ ''^ deTT[niGj,ix)^^^n (28) 



4ei?d Jo 

where (G'jv(x)^^^)^ denotes the Keldysh component of {Gn{x)^^^). In the actual 
calculation, we introduce a parameter 9{x) which is a measure of the proximity effect in DN 
where we denoted ^(0) = 6q in the previous subsections. Using 9{x), Rn{x) can be denoted 
as 

Rn{x) =f3Cose{x) +f2sme{x). (29) 

An{x) and Kn{x) satisfy the following equations, An[x) — — T3^]y(a;)f3, and Kn{x) — 
RN{x)fi{x) — fi{x)A]s[{x) with the distribution function fi{x) which is given by fi{x) — 
foN{x)+f3f3N{x). In the above, fsN^x) is the relevant distribution function which determines 
the conductance of the junction we are now concentrating on. From the retarded or advanced 
component of the Usadel equation, the spatial dependence of 9{x) is determined by the 
following equation 



dx 

while for the Keldysh component we obtain 



D—e{x) + 2ie sm[9{x)] = 0, (30) 



cosh 9imag{x)] = 0. (31) 
14 



with 9imag{x) = Imag[6'(a;)] . At a; = —L, since DN is attached to the normal electrode, 
9{—L)=0 and f3N{~L) = fto is satisfied with 

/to = ^{tanh[(£ + eV)/{2kBT)] - tanh[(£ - eV)/{2kBT)]}, 

where V is the applied bias voltage. Next, we focus on the boundary condition at the DN/S 
interface. Taking the retarded part of Eq. (|26|) . we obtain 

L de{x) . <F> 



Rd dx ^ ^ Rb ^ 
2T„[cos^o(/+ + /-) - sin 60(9+ + g-)] 



(32) 



(2 - r„)(l + g+g. + /+/_) + r„[cos^o(^7+ + g-) + sin^o(/+ + /-)] ' 
On the other hand, from the Keldysh part of Eq. ()26|). we obtain 



-(^)cosh ^0 U^o_= ^ , (33) 



J ^Tn_ go/3jv(0-) 

'° 2 I (2-T„)(l + (?+(7- + /+/-)+T„[cos^o(^7+ + ^?-) + sin^o(/+ + /-)] ? 

Co = T,(1+ I cos^o I' + I sin^o ?) 

X [I (7+ + |2 + I /+ + /_ |2 + I 1 + + g^g_ |2 + I f^g_ _ g^f_ |2] 

+2(2 - TjReal{(l + glg*_ + /;/!)[(cos^o + cos^*)((?+ + g^) + (sin^o + sin^*)(/+ + /_)]} 

+4TJmag(cos^osine*)Imag[(/+ + f-){gX + g*_)]. 
After a simple manipulation, we can obtain /3Ar(0_) 

L J-L COsh^ 0imag{x) 

Since the electric current I^i can be expressed via 60 in the following form 

^ / cosh^[Imag(6'o)](i£, 

we obtain the following final result for the current 



1 r. ft 



e 



Iei = -l de———-^ . (34) 



+ ^ r° 

" <4n> ^ L J- 



dx 

-L COsh^ 9ini.ag(x) 
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Then the total resistance R at zero temperature is given by 

° dx 

2 . ^ - (35) 

L cosh 9imag{x) 

In the following section, we will discuss the normalized conductance arieV) — 
as{eV)/aN{eV) where (75(jv)(eV^) is the voltage-dependent conductance in the superconduct- 
ing (normal) state given by as{eV) — 1/R and aN{eV) — — '^/{Rd + Rb), respectively. 

It should be remarked that in the present circuit theory, Ra/Rb can be varied indepen- 
dently of Tn,i.e. of Z, since we can change the magnitude of the constriction area inde- 
pendently. In the other words, Rd/Rb is no more proportional to Tav{L/l), where Tav is 
the averaged transmissivity and / is the mean free path in the diffusive region, respectively. 
Based on this fact, we can choose Rd/Rb and Z as independent parameters. 

III. RESULTS 

In this subsection, we focus on the line shapes of the conductance where d-wave symmetry 
is chosen as a pairing symmetry of unconventional superconductor. The pair potentials A±. 
are given by Aj. = Aocos[2(0 =]= a)] where a denotes the angle between the normal to the 
interface and the crystal axis of d-wave superconductors and Aq is the maximum amplitude of 
the pair potential. In the above, denotes the injection angle of the quasiparticle measured 
form the x-axis. It is known that quasiparticles with injection angle with 7r/4— | a |<| 
(j) \ < 7r/4:+ I a I feel the MABS at the interface which induces ZBCP. 

A. a = without MABS 

Let us first choose o; = where ZBCP due to the MABS is absent. We choose relatively 
strong barrier Z — 10 (Fig. 1) for various Rd/Rb- For Eth — [see Fig. 1(a)], the 
magnitude of arieV) for | |< Aq increases with the increase of Rd/Rb- First, the fine 
shape of the voltage-dependent conductance remains to be V shaped and only the height 
of the bottom value is enhanced (curve b and c). The V shaped line shape originates from 
the existence of nodes of the ci-wave pair potential. Then, with a further increase of Rd/Rb, 
a rounded bottom structure appears (curves d and e). For Eth = O.OIAq [Fig. 1(b)], the 
magnitude of axieV) has a ZBCP once the magnitude of Rd/Rb deviates slightly from 0. 



„ Rb Rfi 
R^ / 
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The order of the magnitude of the ZBCP width is given by Exh as in the case of s-wave 
junctions js^. When the magnitude of Rd/Rb exceeds unity, the arieV) acquires a zero 
bias conductance dip (ZBCD) (curve e). The quahtative features of hne shapes of (Jr(eV^) 
is different from those in s-wave junctions (see Figs. 1 and 2 in Ref. It should be 

remarked that even in the case of d-wave junctions we can expect ZBCP by CAR as in the 
case of s-wave junction for a = 0. 



a=0 Z=10 
(a) (b) 




FIG. 1: Normalized conductance ax{eV) for Z=10, and a = 0. {a,)Exh = ^o- {]o)ETh = O.OIAq. 
a, Rd/Rb = 0; b, Rd/Rb = 0.1; c, Rd/Rb = 1; d, Rd/Rb = 2; and e, Rd/Rb = 10. 

On the other hand, for much more transparent case with Z = 1 the line shapes of the 
conductance become significantly different. For Exh = Aq [Fig. 2(a)], the magnitude of 
(TT{eV) decreases with the increase of the magnitude of Rd/Rb where the bottom parts of all 
curves are V shaped structures. On the other hand, for E^h = O.OIAq, axieV) has a ZBCD 
for Rd/Rb > 1- Both for Figs. 2(a) and 2(b), the magnitude of arieV) for | el^ |< Aq 
decreases with the increase of Rd/Rb- These features are significantly different from those 
shown in Fig. 1. 

For transparent limit Z = 0, the magnitude of aT{eV) decreases with the increase of 
Rd/Rb (see Fig. 3). For Eth = Aq, arieV) has a broad ZBCP for small Rd/Rb- However, 
with the increase of Rd/Rb, crT{eV) for | eV |< Aq is reduced and becomes nearly constant. 
For Ej^fi = O.OIAq with Rd/Rb = 10, tiny ZBCD appears [curve e of Fig. 3(b)] H- As 
compared to the corresponding case of s-wave junction (see Fig. 4 of Ref. ||38|), ZBCD is 
hard to be visible in (i-wave junctions. 

It is interesting to study how various parameters influence the proximity effect. The 
measure of the proximity effect at the DN/US interface ^q is plotted for Z = and Z = 10 
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with corresponding parameters in Figs. 1 and 3 (see Figs. 4 and 5). For Ra/Rt = 0, 9q = 
is satisfied for any Exh and Z. Besides tliis fact, at £ = 0, 6*0 always becomes a real number. 
These features are consistent with those in s-wave junctions [ssl. First, we study the case 
of Exh/AQ = 1 (Fig. 4) where the same values of Rd/Rb are chosen as in Figs. 1 and 3. The 
real part of is enhanced with an increase in R^/Rb and decreases as a function of e. At 
the same time, the imaginary part of 6q is an increasing function of e for e < Aq. Both real 
and imaginary parts have a sudden change at e = Aq where Aq is the magnitude of the pair 
potential felt by quasiparticles with the perpendicular injection. It is remarkable that the 
magnitude of Real(Imag)(6'o) is reduced with the decrease of Z. Next, we discuss the line 
shapes of 6q for Eth/Aq = 0.01. Real(6'o) has a peak at zero voltage and decreases with the 
increase of e. Imag(^o) increases sharply from and has a peak at about e ~ E^h, except 
for a sufficiently large value of Rd- These features are consistent with s-wave junctions (see 
Fig. 7ofRef. 



38[). Besides this, both real and imaginary parts have a sudden change at 
e = Aq as in the case of Exh/ Aq = 1. Also in this case, the magnitude of Real (Imag) (6*0) is 
reduced with the decrease of Z. This feature can be qualitatively explained as follows. We 
concentrate on the limiting case e = for the simplicity. The magnitude of Oo{0) = 600 is 
determined by the following equation [see eq. 



^ ^< ^(0) >^ COS (PF{<f))d<P 

Rb f:i%cos(f)T{<f))d(j)' 



(36) 



2T„, cos 6*005 



2-Tn + Tn sin 9qo6 ' 



T(0) 



0=0 Z=1 
(a)ETh=Ao (b) Eti^O.OIAq 




-10 1-1 1 
eV/Ao eV/Ao 



FIG. 2: Normalized conductance arieV) for Z=l, and a = 0. (a)ETh = ^o- {^)ETh = O.OIAq. a, 
Rd/Rb = 0; b, Rd/Rb = 0.1; c, Rd/Rb = 1; d, Rd/Rb = 2; and e, Rd/Rb = 10. 
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cx=0 Z=0 



(a) Eth=Ao (b)ETh=0-01Ao 



2 
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CD 
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-1 1-1 1 
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FIG. 3: Normalized conductance aT{eV) for Z = 0, and a = 0. {a)Exh = ^o- (b)-ETh = O.OIAq. 
a, = 0; b, = 0.1; c, Rd/Rb = 1; d, ii^/i?;, = 2; and e, fi^/iib = 10. 

since = = 6 and = = are satisfied with 5 = 1(— 1) for — 7r/4 < (p < tt/4 
(7r/4 <| |< 7r/2). The sign change nature of 6 originates from rf-wave profile of the pair 
potential. This sign change reduces the magnitude of the right hand side of eq. and 
the resulting 6^00 is small. For the case of large magnitude of Z, the degree of the reduction 
due to the sign change of 6, i.e., F{(j)) is not significant. For large magnitude of Z due to the 
existence of the factor T„ (T„ << 1) proportional to cos^0, only the small value of can 
contribute to the integral of numerator where 5=1. This is the reason why the obtained 
measure of 600 for Z = is much smaller than that for Z = 10. 

Although the magnitude of 6q, i.e. the measure of the proximity effect, is enhanced with 
increasing Rd/Rb, its influence on axieV) is different for low and high transparent junctions. 
In the low transparent junctions, the increase in the magnitude of Oq by Rd/Rb can enhance 
the conductance axieV) for eV ~ and produce a ZBCP, whereas in high transparent 
junctions the enhancement of Oq produces the ZBCD. However, the amplitude of dip is 
reduced since the magnitude of 6^00 becomes small due to the sign change of F{(f)). 

B. a / with MABS 

In this subsection, we focus on axieV) and 6*0 for a 7^ (0 < a < 7r/4). First we focus 
on a = vr/S, where MABS is formed for 7r/8 <| |< Svr/S. In the low transparent case, 
i.e. Z=10, axieV) has a ZBCP due to the formation of MABS at the DN/US interface. 
The height of ZBCP is reduced with the increase of Rd/Rb (see Fig. 6). Contrary to the 
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FIG. 4: Real(upper panels) and imaginary parts(lower panels) of are plotted as a function of 
£. Z = 10 (left panels) and Z = (right panels) with Eth/^o = 1 and a = 0. a, Rd/Rb = 0.1; b, 
= 1; and c, Ra/Rb = 10. 
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FIG. 5: Real(uppcr panels) and imaginary parts(lower panels) of Oq are plotted as a function of 
e. Z = 10 (left panels) and Z = (right panels) with Eth/^o = 0.01 and a = 0. a, Rd/Rb = 0.1; 
b, Rd/Rb = 1; and c, Rd/Rb = 10. 

corresponding case of a = (see Fig. 1), arieV) is almost independent of Exh- For Z = 
(see Fig. 7), arieV) has a broad ZBCP both for Eth = Aq and Eth = 0.01 Aq. With the 
increase of Rd/Rb, only the magnitude of axieV) is reduced and ZBCD does not appear. 
It is also interesting to see how Oq is influenced by various parameters. In Fig. 8, line 
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FIG. 6: Normalized conductance arieV) for a = tt/S. {a)ETh = Aq. {h)ETh = O.OIAq. a, 
Rd/Rb = 0; b, Rd/Rb = 1; and c, Ra/Rb = 10. 

a=7i:/8 Z=0 



(a) ETh=Ao (b) ETh=0-01Ao 
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FIG. 7: Normalized conductance arieV) for Z = 0, and a = tt/S. {a)ETh = Aq. {h)ETh = O.OIAq. 
a, Rd/Rb = 0; b, Rd/Rb = 1; and c, Rd/Rb = 10. 

shapes of 9o for Eth/ — 1 are plotted for various parameters. Real(Imag) (^o) has a 
complex e dependence. For Z — 10, the magnitude of Real(^o) is at e ~ is drastically 
suppressed and is an increasing function of e contrary to the case of a 0. Both Real(^o) and 
Imag(^^o) have a peak around e ~ O.TAq ~ Aq cos(2q;) where Aq cos(2«) is the magnitude 
of the pair potential felt by quasiparticles with perpendicular injection. For Z = 0, s 
dependence of Oq is qualitatively similar to the corresponding case of a = (see Fig.4) since 
the role of MABS is not important. For Eth/^o — 0.01, the magnitude of Real(^o) at e ~ 
is suppressed (see Fig. 9) for Z — 10. Both Real(^o) and Imag(^o) have a peak around 
e ~ 0.7Ao as in the case of Eth/^o — 1- On the other hand, for high transparent case, i.e., 
Z = 0, Real(^o) has a peak at e = and decreases with the increase of s. Imag(6'o) increases 
sharply from and has a peak at about e ~ Eth, except for a sufficiently large value of Rd- 
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These features are qualitatively consistent with a — 0. (see Fig. 5). 
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FIG. 8: Real(upper panels) and imaginary parts(lower panels) of 9o are plotted as a function of 
£. Z = 10 (left panels) and Z = (right panels) with Eth/^o = 1 and a = 7r/8. a, Rd/Rb = 0-1; 
b, Rd/Rb = 1; and c, Rd/Rb = 10- 
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FIG. 9: Real(upper panels) and imaginary parts(lower panels) of 9o are plotted as a function of e. 
Z =10 (left panels) and Z = (right panels) with Exh/^Q = 0.01 and a = tt/B. a, Rd/Rb = 0.1; 
b, Rd/Rb = 1; and c, Rd/Rb = 10. 

In order to understand low energy profiles, we focus on the case with e 

determined by the following relation for a ^ 0. 

9oo ^ < F{<P) > 
Rd Rb 



= 0. 000 is 



(37) 
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2-T(</.)+r(</.)sin6»oo w ^| |^ ^ 

-2tan6'oo n/A- a <\ (f) \< n/A + a 

-2T{4,) cos epo 



-^^ wcos!^oo ttM + a <| (A |< 7r/2 

After simple algebra, we obtain 

^00 fsi + /52 + /53 



7r/4— « 



(38) 



/ o /9. • / /■ 2T(0)cosgoo 

/51 = -2V2 tan 6*00 8111 a, 752= / cos0- — , 

Jo 2 - T(0) + T(0) s 



T(0)+T(0) sin ^00 

f 2T(0)cosgoo , r'\(^\ 

fs3 = - COS0- — — . #, /tv = / T COS 

iV4+« 2 - T(0) - T(0) sm ^00 Jo 
In the above, /51 denotes the contribution to < F{(j)) > from the unconventional channel 
with MABS and /s'2 and fs3 denote that from the conventional channel without MABS. For 
low transparent case, T(0) << 1, the magnitude of fsi dominates over those of /s'2 and fs3 
and 6*00 is determined by 

^00 ^ -2^2 tan goo sin a ^^^^ 

Rd RbfN 

Then the resulting ^oo is reduced to be almost zero. As seen from this, MABS and proximity 
effect strongly competes each other. While for high transparent limit, T(0) = 1, the mag- 
nitude of fsi becomes the same order as those of fs2 and fs3- Then from the contribution 
by conventional channel, i.e. fs2 and fss, the magnitude of 6'oo is much larger than that for 
Z = 10. 

For a = 7r/4, where fs2 = fss = is satisfied, only the unconventional channel fsi can 
contribute to < F{(f)) >. N ot on ly 6*00 but 6*0 for any e is exactly zero. Then the total 
resistance R can be given by 

R = +Rd = Rr,=o + Rd (40) 



and the resulting axieV) is given by 

Rd + Rb 



arieV) 



RRa=0 + Rd 

One of the typical examples is plotted in Fig. 10. 

The effect of Rd is significant for the resulting aT{eV). For the actual quantitative 
comparison with tunneling experiments, we must take into account the effect of Rd- 
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FIG. 10: Normalized conductance axieV) for Z = 10, Eth/^o = 1 and a = 7r/4. a, Ra/Rb = 0; b, 
Rd/Rb = 1; and c, Rd/Rb = 10- In this case, arieV) is completely independent of the magnitude 
of Exh- 

C. a dependence of zero- voltage conductance 



Finally, we study the dependence of (7^(0) on the angle a. In this case, (7^(0) is inde- 
pendent of Exh- For all situations shown in Figs. 11 to 13, ^oo = is satisfied for a — n/A, 
due to the complete absence of proximity effect where only the unconventional channel with 

MABS exists. For Z = 0, crr(O) is almost constant with the change of a although ^oo is a 
decreasing function of a (see Fig. 11). 

z=o 
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(b) 
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FIG. 11: Normalized conductance at zero voltage (Tt(0) for various a with Z = 0. In this case 
(Tt(0) is independent of Et^. a, Rd/Rb = 0; b, Rd/Rb = 0-1; c, Rd/Rb = 1; and d, Rd/Rb = 10- 

For Z — 1, (Tt(0) is an increasing function of a, while ^oo is a decreasing function (see Fig. 
12). For Z = 10, (7r(0) is enhanced much more rapidly as compared to the case for Z = 1, 
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while with the increase of Ra/Rb, ct(0) becomes nearly constant (see Fig. 13). As seen 
from Figs. 12 and Figs. 13, the influence of Rd/Rb on the (Jt(0) is significantly important. 
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0.1 0.2 0.1 0.2 

(X/31 a/31 



FIG. 12: Normalized conductance at zero voltage (7^(0) for various a with Z = 1. In this case 
(Tt(0) is independent of Eth- a, Rd/Rb = 0; b, Rd/Rb = 0.1; c, Rd/Rb = 1; and d, Rd/Rb = 10- 
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FIG. 13: Normalized conductance at zero voltage (Tt(0) for various a with Z = 10. In this case 
aT(0) is independent of Eth- a, Rd/Rb = 0; b, Rd/Rb = 0.1; c, Rd/Rb = 1; and d, Rd/Rb = 10. 

In order to understand these features we look at 1^0 in detail. In general for e = 0, < Ibo > 
can be expressed by 

< Ibo >= + t ^ (41) 

-'n 

»7r/4— a 



, . . 2n r r , 2T(c/))[T(0) + {2-T(0)}singoo] 
/m = 2V2smQ;sec ^oo, ^2 = / cos0 — — ... , . „ — ^ , 

, 2T(0)[T(0)-{2-T(0)}singoo] T^^,,, 

^63= / COS(/) — — — . , /„= / r(0)cOS0d0. 

I 2 - T((/.) - r((/)) sm ^00 r Jo 



' 7r/4+a 

Jfei denotes the contribution from the unconventional channel while Ib2 and Jfes denote those 
from the conventional channel. For Z = 0, the integral can be performed analytically. For 
a — < Ibo > becomes 2sec^^oo[l — (v^ — 1) sin^oo] for a = and 2sec^^oo for a — 7r/4. 
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Since the order of ^oo is at most 0.1 as seen from Fig. 11(b), the difference between the two 
cases is smalL Then we can expect that crr(O) is almost constant as a function of a. On the 
other hand, for Z 7^ 0, the contribution from the unconventional channel becomes significant 
and the resulting < /{,o > can be approximated to be 

sin ttsec ^00 /-^n 

It is an increasing function of a, and the resulting o"t(0) is also an increasing function of a 
(see Figs 12 and 13). 



IV. CONCLUSIONS 



In the present paper, detailed theoretical investigation of the voltage-dependent con- 
ductance of diffusive normal metal / unconventional superconductor (DN/US) junctions is 
presented. We have provided the detailed derivation of the expression for the matrix current 
presented in our previous paper |l08| . For the reader's convenience, we explicitly present 
the retarded and the Keldysh parts of the matrix current for the case when the US has a 
singlet parity. Applying these expressions to DN/d-wave junctions, we have obtained the 
following main results. 

1. There are two kinds of ZBCP, i.e. ZBCP due to the coherent Andreev reflection (CAR) 
by proximity effect in DN and that due to the formation of midgap Andreev bound states 
(MABS) at interfaces of (i-wave superconductors. ZBCP frequently appears in the line 
shapes of OTieV). For low transparent junctions with small Thouless energy Et^ we always 
expect ZBCP independent of a except for the large magnitude of Rd/Rb- 

2. The nature of ZBCP due to the MABS and that by CAR is significantly different. The 
corresponding crT'(O) for the former case can take arbitrary values exceeding unity. On the 
other hand, o"t'(0) for the latter case never exceeds unity. The width of the ZBCP in the 
former case is determined by the transparency of the junction while the width for the latter 
case is determined by the Thouless energy. These two ZBCPs compete each other sinc e the 
proximity effect and the existence of MABS are incompatible in singlet junctions 108 ]. 

3. For the extreme case, a = tt/A, where the proximity effect is absent and the CAR is 
cancelled, axieV) is given by axieV) = {Rb + Rd)l{RR^=o + Rd) with the resistance at the 
interface Rb- 4. Only when a = 0, MABS is absent for Rj, = 0. Then CAR influences 
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significantly on aT{eV), similarly to the case of s-wave junction. When the transparency 
of the junction is sufficiently low, aT{eV) for | eV |< Aq is enhanced with the increase of 
Rd due to the enhancement of the proximity effect assisted by CAR. The ZBCP becomes 
prominent for Eth << Aq and Rd/Rt < In such a case, with a further increase of Ra/Rb 
the ZBCP changes into a zero bias conductance dip (ZBCD). 

5. We have clarified various line shapes of the conductance including ZBCP. The obtained 
results serve as an important guide to analyze the actual experimental data of the tunneling 
spectra of high Tc cuprate junctions. We want to stress that the height of ZBCP is strongly 
suppressed by the existence of DN and the resulting crr(O) is not so high as obtained in the 
ballistic regime j4o|- In the actual fit of the experimental data, we strongly hope to take 
into account the effect of Rd- When the transparency of the junction is low and a 7^ 0, 
the contribution of unconventional channel becomes important and that from conventional 
channel is negligible. In such a case without solving Usadel equation aT{eV) can be simply 
approximated by 

^T{eV) = - -— , Rj,^^, = ——— (42) 

rCR^=o + rid < ^bo > 



< ho >- 



cos ^T{(j))d(j)' 



T(0){1 + T(0) 


r+ 


' - 1] 


r+r_ 






1 + [T(0) - i]r+r_ 


2 



A^ _ A_ h 

Rb 



e+v/i^^A^' " 2e2jV2^rf0T(0)cos0' 
with A-i- = Aq cos[2(0=i=a)] and e = eV. This expression is a convenient one for the fit of the 
experimental data. However, for the quantitative discussions including much more general 
cases, one must solve Usadel equation as was done in the present paper. It is an interesting 
future problem to compare the present results with experiments since recent experi ment al 



results show the existence of the mesoscopic coherence in high Tc cuprate junctions [11 

There are several problems which are not discussed in the present paper. In the present 
study, we have focused on N/S junctions. The extension of the circuit theory to long 
diffusive S/N/S junctions has been performed by Bezuglyi et al 3^. In S/N/S junctions. 
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e Andreev ref l ections pr oduces the subharmonic gap structures 
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complex as compared to N/S junctions. Moreover, in S/N/S junctions with unconventional 
superconductors, MABS leads to the anomalous current-phase relation and temperature 
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dependence of the Josephson current [SOj. An interesting problem is an extension of the 
circuit theory to S/N/S junctions with unconventional superconductors. 

There are two kinds of ZBCPs considered in the present paper. We expect that the 
response to the magnetic field should be significantly different in these two cases. The 
ZBCP originating from MABS is rather robust against the magnetic field while that from 
CAR is much more sensitive. We want to clarify this feature in actual calculations. 

In the present paper, since we follow the quasiclassical Green's function formalism, the 
impurity scattering is taken into account within the self-consistent Born approximation. It 
is a challenging problem to study the weak localization effects. 

The authors appreciate useful and fruitful discussions with J. Inoue, H. Itoh, Y. Asano 
and I. Shigeta. This work was supported by the Core Research for Evolutional Science 
and Technology (CREST) of the Japan Science and Technology Corporation (JST). The 
computational aspect of this work has been performed at the facilities of the Supercomputer 
Center, Institute for Solid State Physics, University of Tokyo and the Computer Center. 

V. APPENDIX 

The matrix current is expressed as 

/ = ^Tr„„,[S^^i] = ^Tr„„,[S^^2]. (43) 

To find the matrix current, we have to evaluate ^1(2)- For this purpose, we shall consider the 
behavior of G'^'{x,x') in the isotropization zone in DN side {—Li < x, x' < —L2) and in 
the ballistic zone of right side (x,x' > 0). In the isotropic zone, since the effect of impurity 
scattering is dominant and H can be approximated to be —T,imp{x) = Gi/{2Timp) for x < 0. 
Gi is the Keldysh-Nambu Green's function in DN at x = —Li with ^1 >> Li >> VTi^p and 
^1 >> L1—L2 >> VTimp where ,^1 = ^J D/2tiT is the coherence length of the Green's function 
in DN. Due to this condition, Gi can be approximated to be Gi = G^i—Li) ~ G^^O^), 
where G^ix) obeys the Usadel equation in DN. The Green's function (^^""'(x, x') is expressed 
by 

G--'(a;,a;') = P{x){gi + sign(x - x')S")P(-x') 
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P{x) = -^={eMx/{2vnrimp)]{l - S^Gi) + exp[-x/{2vnrimp)]{l + S^Gi)} (44) 

with 



, \ - / i , 

GiJ yo -1 

To ensure that G^"''{x,x') does not grow with decreasing x and x' in the isotropic zone in 
DN, we shall require 

(S^ + Gi)(S^-^i) = 0, (46) 

(E^ + ^i)(E^-G'i) = 0. (47) 

On the other hand, in the US side {x > 0) it is not simple to obtain G^^' {x, x') since it has 
directional dependence. G"^' {x,x') is given by 

Gl^'ix.x') = P{x){g2 + sign(x - x')Y.')P{-x'), 



P{x) = Pi{x) + P2{x), (48) 
A(2)+ 6 



Pl(2)- 



5 , , 1 7±(1t4±) /o[7±(lT4±)-7i(lT^±)] 
-'lit l,-^ j — — I I 

\^ 7^(lip^2±) 

^ , , I ( 7±(1±4±) /o[7±(l±4±)-7i(l±^±)] 

2y/2vni \ 7^(1 ±i2±) 



with 



and 



cxp [—i ' ^ x] 



7± = exp z x\. 

hVFx 
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In the above R2±, Ai± are retarded and advanced component of Keldysh-Nambu Green's 
function G2± at the interface of US where ± denotes the direction of motion along x- 
axis. 6*2+ and G2- are given by G2± = Gs±i^+) where Gs±{x) is a quasi classical Green's 
function in US. It does depend on the direction of motion a. Here we neglect the spatial 
dependence for simplicity and we assume G2± = Gs±{x) = Gs±{oo). Since ■y± and ■y'^ are 
glowing functions with the increase of x, the term, which includes this component should 
be eliminated by multiplying (^2 + sign(a; — For convenience, we denote 





g2 + sign(x — x )S 

In order to eliminate the divergence terms with the increase of x{x > 0) 
7+(l - R2+) /o[7+(l - 4+) - 7;(1 - ^2+)] W « c \ ^ ^ ^ 

7;(i-^+) ) V V 

should be satisfied. Then following four equations must be satisfied for any x 

7+(l - R2+)a + /o[7+(l - R2+) - 7:(1 - = 0, 

7;(l-i2+)& = 0, 

7;(l-i2+)c^ = 0, 
7+(l - R2+)d + /o[7+(l - R2+) -llil- A2+)]d = 0. 

Thus we obtain 

(1-^+) ) \b d) 

Applying the similar discussions for other components, we also obtain the following equa- 
tions, 

^ E, \ f A C 



0, 



\ B D 



with 



'l-R2±) fo[-R2±+A2±)] 

a-A2±) 



The resulting equation is identical to eq.(30) of Ref. (37 1. 



(S^-G'2)(^2 + S^) = 0. (49) 
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Applying similar discussion for the case with increasing x' [x' > 0), we can also obtain 



(^2-S^)(S^ + G'2) = 0. 

To find explicit expression of (72, we multiply Eq. EHlby M^" from the left and by M from 
the right. By introducing the matrix Q = M^M, we derive 

92 = {QG2 + Gi)-^{2Q + (Gi - (50) 

with 



Go 



G2+ 
G2- 

To evaluate matrix current /, we rewrite 92 in the basis composed of eigenvectors of Q and 
Keldysh-Nambu indices. For each eigenvector c„ with eigenvalue g„ > 1, the vector S^c„ is 
also an eigenvector of Q with eigenvalue 

The resulting ^2 is given by 

qnH+ + G, qnH^ \ V Qni2 - H^) Gi-qnH+ ^^^^ 

with H± = (G2+ ± G2-)/2. Then the matrix current / can be expressed as 

n 



Lo = [H^ - {qnH+ + Gi)Hz\H+/qn + Gi)]-\2 - H.) 

+ [H- - {H+/qn + Gi)Hz\qnH+ + Gi)]-\2 - H^) 

+ [{H+/qn + Gi) - H-{qnH+ + Gi)-'H-]-\Gi - H+fq^) 

+ [{qnH+ + Gi) - H.{H+/qn + Gi)-^i^_]-^(Gi - qnH+) (53) 
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Applying the following identity, 

Hl + Hl = l H+H- + H-H+ = 
J„o can be written as follows, 

J„o = {[Gi, HZ^]qn - HZ'H+ + G^Hz'H+Grql)-\Gr - Hz'G^)qn + G^HZ^H^G^ql] 

+ ([61, HZ^yqn - HZ^H+lql + G^HZ^ H+Gi)-\Gi - G,Hz')/qn + Hz'H+/ql] 
+ ([61, HZ']/qn - HZ'H+ + GiHz'H+GJql)-'[iGi - HZ'G,)/qn + GiHz'H+G,/ql] 

+ Hz']qn - Hz'H+ql + G,Hz'H+G,)-'[{G, - G,Hz')qn + Hz'H+ql] (54) 
The eigenvalue g„ and the transparency of the junction T„ satisfy the following relations, 
4gn/(l + qn)' = r„, (1 + g2)/(l + g„)2 = (2 _ T„)/2 

We introduce Tin 

T 

In 



2 - T„ + 2v/r^7; 

and the resulting I^o becomes 

Ino = -D-'[Tin{2Gi - [HZ\ GiU) + Hz'H+ + TIG^HZ^H+Gi] 

+GiD-'Gi[Tini2Gi - [HZ\ Gi]+) + KHz'H+ + G^HZ^H+Gi] (55) 

with 

D = —Tin[Gi, HZ^] + HZ^H+ — T^^GiHZ^H^Gi 
Ino is also represented as follows liosl ]. 

InO = 2[G'l, -Bn] 

with 

Bn = {-Tin[Gi,HZ^]+HZ^H+-TlfiiHZ^H+Gi)-\Tin{l-HZ^ (56) 
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